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1 Introduction 



Analytic solutions to the Navier-Stokes equations are few and far between. In [T] a new 
solution to the relativistic, conformally invariant Navier-Stokes equation was constructed. 
The solution is a generalization of Bjorken flow [2] and describes a boost-invariant medium 
expanding both longitudinally and radially, which makes it an attractive toy model for heavy 
ion collisions. We will refer to this solution as the 5*0(3) x 5*0(1, 1) x Z2-invariant solution, 
or the 50 (3)-invariant flow for short. The former title describes the symmetries used to 
construct the solution: like Bjorken flow, it is boost invariant, and it has a Z2 symmetry 
generated by reflections through the collision plane; but instead of being translationally 
invariant in the transverse directions it has an 5*0(3) conformal symmetry. The explicit 
form of the inviscid 5*0(3)-invariant flow, and a simple derivation of it, can be found in 
section |2l 

Conformal invariance was a crucial assumption in the construction of the 5*0(3)-invariant 
flow. Conformal invariance implies that the dynamics is invariant under a Weyl rescaling of 
the metric: 

il^'^g^^ , (1) 

where Q is allowed to vary across spacetime. The Minkowski metric, the de-Sitter (dS) 
metric, and the anti-de Sitter (AdS) metric are all related to one another through Weyl 
rescalings. A Weyl rescaling locally preserves angles but not sizes. When considering the 
hydrodynamic evolution of the presumed quark-gluon-plasma (QGP) at RHIC and the LHC, 
it is a reasonable approximation in the flrst few fm/c to set the pressure p equal to e/3, 
and to set bulk viscosity equal to 0. These are precisely the conditions guaranteed by the 
tracelessness of the stress tensor, which is the hallmark of conformal invariance. In other 
words, p = e/3 and C = means that the symmetries of relativistic hydrodynamics are 
enlarged to include conformal symmetry. We discuss conformal hydrodynamics at greater 
length in section |3| 

Conformal symmetry turns out to be a powerful tool for generating non-trivial solutions 
to the Navier-Stokes equations, both exact and approximate. Our derivation in section |2] 
of the inviscid 5*0 (3) -invariant flow starts with a stationary fluid in the geometry dS^ x R. 
Owing to the simplicity of the solution in dSs x R coordinates, in section |4| we are able to 
generalize it to arbitrary space-time dimensions and to include the effects of a non-vanishing 
chemical potential and second order viscous corrections. In section [5] we explain how to 
rewrite this solution as a dual time-dependent AdS black hole geometry. 
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In section [6] we derive the equations that govern arbitrary hnear perturbations of the 
5'0(3)-invariant flow. A key question, which we study in section [7| is whether there are 
instabihties. The answer is that instabihties do exist, but for choices of parameters chosen 
to match approximately with RHIC collisions, the instabilities occur at times earlier than the 
thermalization time. Because of the recent interest in higher order moments of the flow of real 
quark-gluon plasmas |3l IH El El |7l [Hj , we also consider, in sections [T] and [sj explicit solutions 
of the equations governing linear perturbations of the S'0(3)-invariant flow. These solutions 
are classifled by their transformation properties under 5*0 (3). We also explain in section |8] 
how perturbations to the initial state can also be classified in terms of representations of 
S0{3). 

2 A simple derivation of the 50 (3) -invariant flow 

Inviscid, conformal, relativistic hydrodynamics in 3 + 1 dimensions is based on the following 
form for the stress tensor: 

T^u = (^U,,Uu + ^Pfiu (2) 

where e is the energy density and 

Pfiu = u^Uu + g^u (3) 

is the tensor that projects onto the local rest frame of a fluid element. The 4-velocity is 
subject to the constraint 

^^X^- = -l- (4) 

The relativistic Euler equations are 

= 0. (5) 

Consider coordinates (r, r], 0) for R'^'-'^, such that 

ds"^ = -dr^ + r'^drf + dx\ + x\d(t)^ . (6) 

This coordinate system only covers a part of R^'^, which is the region in the causal future 
of the collision plane at r = 0. We will call this region the future wedge. A solution to 
equations ([s]) was studied in ^ which is deflned on the future wedge of R^'^ and takes the 
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following form: 



Ur = — cosh k{t,x±) m_l = sinli K(r, = = (7) 

where the transverse or radial velocity is given by 

v± = tanh k{t, x±} = = , (8) 

Ur 1 + q^T^ + 

with g an arbitrary dimensionful constant with units of inverse length, and 

60 (2g)«/^ 



[1 + 2g2(r2 + xl) + q\T^ - xlff^ 



(9) 



with eo a dimensionless integration constant. The solution ([7])-([9]) was obtained almost 
entirely through symmetry considerations. However, the methods used were less than trans- 
parent because the symmetries in question are not isometrics of R'^'^, but instead conformal 
isometrics. The purpose of this section is to re-derive this solution using methods that make 
the symmetries more manifest. A more detailed exposition of our new method, including 
viscous corrections, can be found in sections [3] and |4} 

The key step is to make a coordinate transformation combined with a Weyl rescaling of 
the metric which promotes the S'0(3) conformal isometry to a manifest isometry. Explicitly, 
the Weyl rescaling is given by 

where ds"^ is the line element The coordinate transformation we use is defined through 
the relations 

. + 2qx± 

smhp = tane = 11 

zqr 1 + q^T^ — q^Xj_ 

so that the rescaled line element takes the final form 

ds^ = -dp^ + cosh^ p {de^ + sin^ Odcp'^) + dif . (12) 
We emphasize that the metric ds^ is not the standard fiat metric on R^'^ because of the 



overall factor of 1/r relative to ds . As it turns out, ds is the metric of dS^ x R. In (12) 



we have passed to standard global coordinates on dS^^. While in (11) is the same angle 



as in ^ (i.e. the angle around the beam-line), 6 in (11) has nothing to do with the polar 
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angle {} = sin~^tanhr7 which vanishes at mid-rapidity; instead, as (11) shows, {p,0) provide 
an alternative parametrization of the directions {t,x±). 

The (p, 6, (f), rj) coordinate system have symmetries that were not manifest before: namely, 
the rotations of the sphere parametrized by {9, (p). These symmetries comprise the promised 
5*0 (3). Its generators are isometries of dS^ x R, but (aside from rotations in the (p direction) 
only conformal isometries of R^'^. It is obvious how to construct an S'0(3)-invariant velocity 
profile in the {p,9,(j),ri) coordinates: 

Up = -1 ue = u^ = Un = 0. (13) 



Clearly, there is no reference to any specific dynamics in (13): we are just describing a fluid 
which is at rest in the dS^ x R geometry. To get back to a velocity profile in the future 
wedge of R^'\ we need to go back from the (p, 9, 0, rj) coordinate system to the (r, x±, (p, rj) 
coordinate system and carry out an appropriate Weyl rescaling of the velocity field: 

% = -|^^- (14) 
where = [p, 9, (p, rj) and x^ = (r, t], x±,(p). The explicit factor of r on the right hand side 



of (14) follows from the unit norm constraint, 

u^^u^g^, = -1 = uVg^, . (15) 



Now we can recover dSl) by plugging (13) into (14): 



(16) 



The partial derivatives in (16) are to be computed starting from the defining equations (11). 
Note that, in common with Bjorken flow, it is assumed that = = 0, i.e. there is 
rotational invariance around the beam-line and also boost invariance. 

We can go further and obtain the energy density of the 5'0(3)-invariant flow as follows. 
In the ^5*3 x R frame, the energy density e must be a function only of de Sitter time p. In 
the absence of a chemical potential the equation of state implies that the entropy density, 
s, is proportional to e"/^. At some fixed de Sitter time p, the total entropy per unit rapidity 
is given by (47r cosh^ p)s: this product is the entropy density times the volume of the S'^ 
parameterized by {9,(p). Because there is no viscosity, the total entropy can't change with 
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p. Thus, we conclude that 

:^(e=^/^coshV) =0, (17) 
dp ^ ' 



whose solution is 



e = eo(coshp) 



where eo is a constant. We use the transformation law 

(19) 



to recover the energy density e on the future wedge of flat R ' . To understand the rescaling 



(19), recall that lengths and times in R^'^ are scaled by r relative to dS^ x R. Volumes 
in R^'^ should therefore be scaled by relative to corresponding volumes in dS?^ x R, and 
energies should be scaled by 1/r; thus energy density should be scaled by 1/t^, as we have 



done in (19). We give a more formal derivation of ( |19[ ) in section |3j 

In summary, the seemingly complicated solution ((7])-([9]) can be recovered through a con- 
formal transformation from isentropic evolution of a static fluid in ^5*3 x R. The simplicity 
of the ^5*3 X R picture invites generalizations, and we will devote the rest of this paper to 
exploring some of them. Shear viscosity (already treated in p]), higher-derivative corrections 
to hydrodynamics, finite chemical potentials, and anisotropics are all susceptible to partially 
analytic treatment in this framework, as are generalizations to other dimensions. 

A peculiar feature of the dS^ x R conformal frame should be mentioned before going 
further: the coordinates (r, a;_L,(/)) cover only half of the full ^5*3 geometry, whereas the 
coordinates (p, 6', 0) cover the whole. Moreover, the half covered by the coordinates (r, 0) 
is a contracting geometry. What this means is that the distance between two points at fixed 
coordinates in the transverse plane decreases as r increases. Explicitly, this distance is ^ if 
X is the coordinate distance in the transverse plane. In appendix |A] we provide more detail 



about different ways of putting coordinates on de Sitter space, and in figure 10 we show how 
surfaces of constant r cut across de Sitter space. 



3 Conformal transformations and hydrodynamics 

Having presented in section [2] a simple example illustrating our new method of finding 
solutions to relativistic conformal hydrodynamics, we aim in this section to give a more 
systematic summary of the conformally invariant Navier-Stokes equations, in preparation 
for exploration of more complicated fluid flows in later sections. The material in this section 
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is not new. Much it is part of the classical literature on hydrodynamics: see for example [U]. 
More recent literature includes [ini [HI [ISl fT3] . 

To fix notation, let's re-express the Weyl rescaling ([T]) as a relation between two metrics 
on the same ci-dimensional spacetime: 

ds' = Q^ds^ . (20) 

Recall that the conformal factor Q is allowed to depend on the spacetime coordinates, but 
should be everywhere non-zero and non-singular. Note that a Weyl rescaling is not a co- 
ordinate transformation: Indeed, we can use the same coordinates for the metrics ds"^ and 
ds^. 

It will be important for us to track how various quantities transform under Weyl rescal- 
ings. In general, we say that a quantity X transforms homogeneously with weight ax under 
a Weyl transformation if 

X = Q-"^X. (21) 
X may be a scalar quantity or a tensor, and we will refer to ax as its conformal weight. We 



will abbreviate (21) as 

[X] = ax . (22) 
Note that the conformal weight ax depends on the index structure of X: For example, we 



read off from (20) that [gfj,u] = —2, but the inverse metric has weight [g'^'^] = +2. 

Now let's discuss the conformal properties of the energy momentum tensor and a con- 
served charge current. The energy momentum tensor is given by 

2 6\nZ , , 

(T^'^) = (23) 

V-fl' ^9iLv 

where Z is the partition function of the theory. In a d dimensional conformal theory, and in 
the absence of conformal anomalies]^ Z is invariant under conformal transformations. Thus, 

[T^n = d^2. (24) 



One can check that ( 24 ) together with tracelessness of the stress tensor implies that energy- 



momentum conservation is a conformal-frame independent statement. A similar relation can 



^In practice, quantum conformal theories are anomalous |14j . However, when considering hydrodynamics 
as a derivative expansion this anomaly appears only at order d (for even d). Thus, to the order we are 
considering in this work the conformal anomaly can be neglected. See |10j for details. 
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be established for a conserved current J^: If we want the conservation equation 



= (25) 

to be conformal-frame independent, then we need 

[r] = d. (26) 

In the hydrodynamic approximation, the degrees of freedom of the theory reduce to the 
energy density e, a velocity field normalized such that u^m^ = — 1, and charge densities 
Hi (which one can substitute for chemical potentials, depending on the choice of ensemble)]^ 
For the sake of clarity we will consider turning on only one type of charge, n. The index 
runs from to d — 1. Since g^uU^W^ = — 1 we find that [u'^] = 1. Further, in a frame where 
u^T^" = —euy we can write 

Ti"' = eu^^u" + Tl" (27) 

with u^T^ = 0. Hence, [e] = d. Similarly one can show that [n] = d — 1. Consequently 
[T] = [fi] = 1 where T is the temperature and fi is the chemical potential. The conformal 
weights of quantities to be encountered in the rest of this paper are summarized table [T} 

The transverse component of the energy momentum tensor can be expanded in a gradient 
expansion. For a conformal theory one finds that to first order in gradients, 

T"" = euV + j^P^'' - W'"" (28) 
where the projection matrix P^^ was defined in ([3]) and 



0" 



2V(^n,) . (29) 



Here, brackets denote a symmetric traceless projection onto the space orthogonal to u^. 
More explicitly, 

= \p,aPup (A"^ + A^") - -^P.^P'^^A^p . (30) 
Similarly, in the presence of a single conserved charge n, and in the absence of anomalies. 



^In the case of a spontaneously broken symmetry the gradient of the Goldstone mode is also a hydrody- 
namic degree of freedom and can be interpreted as a superfluid velocity. 
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quantity 


conformal weight 


description 


extra relations 




-2 


metric 


ds"^ = Qfj^ydx^dx^ 




-1 


velocity 


g^^u^Uy = -1 


p 


-2 


nroipcts onto 
local rest frame 




T 


d-2 


stress-energy tensor 


rpH ^ Q V^T,,^ = 0, 

, , 2 5\nZ 

{Tn = — 

V-g dg^,u 


e 


d 


energy density 


T^,u = eui^uy + pP/^y + ... 


P 


d 


pressure 


€ 


T 


1 


temperature 


e oc if - 




d-2 


entropy current 


> 0, S/, = SM^ + . . . 


J, 


d-2 


conserved current 




n 


d-1 


charge density 


Jfj. = nu^ + ... 




1 


chemical potential 


e + p = Ts + iin 



Table 1: Conformal weights and additional relations for tensors of interest. The additional 
relations are valid provided the theory is free of anomalies and there are no Goldstone 
modes. Omitted terms, denoted are higher derivative corrections. Lorentz indices are 
raised using g'^'^, which changes the conformal weight by +2. 
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the associated conserved current J is given by 



(31) 



The relativistic version of the Navier-Stokes equation are simply energy-momentum conser- 
vation: 



4 Solutions with SO{d — 1) x 50(1, 1) x Z2 symmetry 



Consider a conformal theory on R'^"^'^. The Minkowski metric r] has d{d— l)/2 Kilhng vec- 
tors and d + 1 additional conformal Killing vectors. As a solution generating technique, an 
obvious strategy is to impose symmetry under a subset of these isometrics and / or conformal 
isometrics. In order to reduce the Navier-Stokes equations from partial differential equations 
(PDEs) with d independent variables (namely the coordinates on R''"^'^) to ordinary dif- 
ferential equations (ODEs), one needs to impose d — 1 independent symmetry constraints. 
This is exactly what is done in Bjorken flow: starting with R^'^, one needs 3 symmetry con- 
straints to get ODEs, and the independent symmetries used are boost invariance, rotational 
invariance around the beamline, and translation in one of the two directions transverse to 
the beamline. Invariance under translations in the second transverse dimension follows, be- 
cause the generator of such translations is the commutator of rotations and translations in 
the first transverse direction. Similarly, the flow explained in section [2] has three indepen- 
dent continuous symmetries: one for 5*0(1, 1) boost invariance plus two for S0{3) conformal 
symmetry, with symmetry under the third generator of SO {3) following from commutations 
of the other two generators. 

The plan of this section is to seek elementary generalizations of the S'0(3)-invariant flow 
explained in section [2] to arbitrary dimension and non-zero shear viscosity. We will use the 
same trick of passing via a Weyl rescaling from the future wedge of R"^^^'^ to dSd-i x R, 
and then (for the most part) consider fluid flows that are stationary in dSd-i x R. By 
assumption, our flows possess SO{d — 1) symmetry (which acts on the S"*^"^ spatial slices of 
dSd-i), S'0(1, 1) boost symmetry (which acts additively on the R factor of the geometry, 
parametrized by rapidity rj), and a Z2 symmetry under rj — )■ —7]. The way the counting of 
the previous paragraph works out is that SO{d — 1) symmetry can be expressed as — 2 




(32) 



In the presence of charges one must also impose (25). 
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independent constraints, while S0{1, 1) imposes one more independent constraint; thus we 
have precisely the d—1 independent constraints we need to reduce the d-dimensional Navier- 
Stokes equations to PDEs. 



4.1 Switching from ^'^ to dSd^i x R 

Let's start by passing from Cartesian to Bjorken coordinates on R'^"^: 

ds'^ = -dt^ + dz^ + (ix^ 



-dr"^ + r'^drf + dx^ 



(33) 



± 5 



where ^) parametrizes the transverse plane R"^ ^, and r and rj are 

related to t and 2; through 

r = ^ — t'^ tanh rj = - . (34) 

The region r > covers the future wedge t > In this coordinate system the 5*0(1, 1) 
symmetry manifests itself as translations in rj. 

We want to switch to the dS^-i x R conformal frame, but first we want to provide a bit 
more detail about dSd-i- Its line element is 

ds"^ = ^ {-dr^ + dxl) . (35) 

It is helpful to recall that dSd-i can be realized as the locus of points in R'^"^'^ satisfying 
the equation 

d-l 

-(X°)2 + ^(X'")2 = 1. (36) 



m=l 



If we parameterize this surface by 



2qT T 2qT ' 



where q is an arbitrary parameter, then (35) can be recovered as the natural metric on the 
locus (36 ) inherited from the standard flat metric on R*^"^'^. On the other hand, if we express 



an arbitrary point on S unit vector in R<^-\ then we may relate 

X° = sinh p X'" = cosh p r"" . (38) 
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with a resulting line element for dS^-i'- 



ds^ = —dp^ + cosh^ p dil'^__2 



(39) 



where dfl'^_2 is the metric of the d — 2-dimensional unit sphere. Note that p is a time 
coordinate: it is global time in ciS'^-ij^ 

Evidently, the future wedge of R'^"^'^ is conformally equivalent to dSd-i x R: 



ds^ = 



-dr"^ + dx^. 



+ d7f 



[—dp^ + cosh^ p (ifi^_2 + dif\ 



(40) 



We note in passing that the future wedge of R*^ ^'^ is also conformally equivalent to AdS2 x 
rjnj^^ Minkowski metric on R*^"^'^ may be expressed as 



ds^ = cosh^ p 



-dp^ + drf 
cosh^ p 



(41) 



In order to present explicit formulas for how hydrodynamic quantities transform in passing 
from R'^"^'^ to dS^-i x R, it is useful first to parametrize the transverse plane R*^"^ in polar 
coordinates, so that the standard metric on the future wedge of R^'^ can be written as 



ds' 



-dr"^ + r'^drf + dx]_ + x]_dVL 



d-3 5 



(42) 



where dQ'^_^ is the metric on the unit sphere S'^~^, parametrized (in a manner we need 
not specify precisely for present purposes) by d — 3 angles (pi. From here on, our preferred 
coordinates for R'^"^'^ will be 



= (r,X_L,0i, . . . ,0d_3,?7) . 



(43) 



On the other hand, we may express the metric on dSd-i x R as 



df = -dp^ + cosh^ p {dO'^ + sin^ 9 dnl_^) + dr]'^ 



(44) 



and our preferred coordinates on dSd-i x R will be 



= {p,9,(f)i, . . . ,0d„3,r7) . 



(45) 



^In comparing to [T] , it is useful to note that g = — sinh p. 
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The map we explained from the future wedge of R'^ ^'^ to dSd-i x R in equations (|33|) 



(40) leaves the (f)i coordinates alone. The map sends {t,x±) — )■ {p,0) through precisely the 



formulas (11) introduced in section 2 The energy density, charge density, and velocity fields 



transform as follows: 



e 


= r e 


n 




Ur 




u± 


-t( u 
\dx± 




= ^U^i 




= TUr, 



de 

-ue 



as ^ 



dx 



where we consistently use hats to denote quantities in the dSd-i x R conformal frame, while 
unhatted quantities are for R'^~^'^. 

In the rest of this paper, we will work almost exclusively in the dSd-i x R conformal 



frame, as it is understood from (46) how to go back to R ' quantities. 



4.2 The inviscid case 

Having chosen the dSd-i x R conformal frame, parametrized by variables x^ = (p, 6, 
01, . . . , 0d_3, 77), let us now consider the velocity profiles it^ permitted by SO{d — 1) x 
S0{1, 1) X Z2 symmetry. The S0{1, 1) x Z2 symmetry implies = and prevents the 
other components of from depending on rj. For d > 3, the SO{d — 1) dependence, on top 
of 5*0(1, 1) X Z2, leaves only one possible velocity profile: 

Up = -1 (47) 

with all other components set to 0|^ When d = 3 the ue component of the velocity profile does 
not necessarily vanish. We will treat the d = 3 case with non-vanishing angular momentum 



in the 9 direction separately in section 4.6 



Before working out the energy density, let's consider how charge density evolves. Without 



appealing to the explicit form (31) of the conserved charge current J'^, we note that the 



The sign in (47) is fixed by noting from the first line of ( 11 ) that r — > corresponds to p — > —00. Thus 



the future direction is toward more positive r and toward more positive p. 
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SO{d — 1) X S0{1, 1) X Z2 forces all components of to vanish except 



jp = -u jf^^n. (48) 



Plugging (48) into the conservation equation (25) leads immediately to 

d 



whose general solution is 



[(cosh p)'^-2n] =0, (49) 



n = no(coshp)-('^-2) ^50^ 



where no is a constant. Note that the result (50) has no dependence on the diffusion constant 



K, which makes sense because no net diffusion occurs during the flow. 

So far we have made no appeal to dynamical equations of motion: in particular, our 



choice of velocity field (47) should be valid equally for the inviscid and viscous cases, and 
the result (50) requires only symmetry considerations plus the conservation law for charge. 
Therefore, all the analysis so far is valid equally in the inviscid and viscous cases. 

In order to determine the energy density for the inviscid case, we can use the same style of 
argument as in section [2] Entropy is not produced in the absence of viscosity, so the entropy 
density s in dSd-i x R behaves just like the conserved charge h: that is, s oc (cosh p)~^'^~'^\ 
Conformal invariance requires that the energy density e in dSd-i x R obeys an equation of 
state of the form 

e = nd^e{a) (51) 



where 

s 



a = -. (52) 
n 



Conformal invariance tells us nothing about the scaling function e(cr). The crucial point is 
that a is constant throughout the flow, simply because n and s have the same dependence 
on p (up to overall constants) in the absence of viscosity. Thus e oc and combining this 



result with (50) we find 



eo(coshp) , (53) 



where eo is a constant. 
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4.3 Viscous corrections 



In section |4.2[ we saw that the solution of conformal relativistic hydrodynamics respecting 
SO{d — 1) X S0{1, 1) X Z2 symmetry is essentially trivial, even for the most general equation 



of state (51) allowed by conformal invariance in the presence of non-zero charge density. 
The reason for this triviality is that the evolution of entropy density in de Sitter time p 
exactly tracks the evolution of the charge density. Viscosity will spoil this. In order to make 
progress, we need a clever parametrization of shear viscosity: 

fj = hl\a)fi{a) , (54) 



where fj is another dimensionless function of a. Because fj is arbitrary, the form (54) is the 
most general one allowed by conformal invariance. It seems that it would be simpler to omit 
the factor e'(cr), but including it is convenient because then e((j) factors out completely from 
the hydrodynamic equations for the stationary fluid in dSd^i x R. These equations boil 
down to 

1 da 2{d-2)^_^. 



, n d-i tanh p, (55) 
rj{a) dp d — 1 

which is separable. A slight additional generalization is helpful: we need not insist that a is 
precisely the ratio s/n; instead, it can be some function of this ratio without changing any 



aspect of the calculations we have done, provided we employ the parametrizations (51) and 



(154|). 

To proceed further, let us define 



^ / X r , / , 1 2 sinh> ^ /3 2c; - 1 5 , , \ 

Fd{p) = J ^'^ (coshr)d-i tanh r = — - — 2F1 (^-, ^ ; -; -smh p 1 . (56) 

For any fixed d, the map p — ?■ Fd{p) is a smooth bijection of R to R, and its slope is positive 



except at p = 0. Moreover, it is odd under p — —p. The general solution of (55) can be 
readily expressed in terms of F^: 

-'^'-'^K^'FAp), (57) 



rj{u) d — 1 
where ctq is the value of a at de Sitter time p = 0. 



The solution (57) is slightly abstract. Let us work out an example to show its utility. 



Consider the case where neither the energy density nor the viscosity depend significantly on 
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the conserved charge. Then we define 



and fj = HqT 



d-l 



(58) 



where T is proportional to the temperature, and Hq is a dimensionless constant. Hq must re- 
main arbitrary if we want to capture the dynamics of arbitrarily viscous fluids. The conserved 



current is now treated in a probe approximation: (50) still holds, but by assumption, n does 



not enter into the stress tensor in any way. We may nevertheless employ the parametrizations 



(51) and (54), for example with 



e~(a) = a° 



Ho 
d 



T 



a 



(59) 



Then it is easy to see that (57) becomes 



f{p) = (coshp) 



2{d-2) ■ 



(60) 



where Tq is the de Sitter temperature at de Sitter time p = 0. Note that n cancels out of the 
final result (60), as it must. Using (46) we can map (60) to the energy density on R'^"^'^. 
We find 



(61) 



For d = 4, the solution (61) reduces to the one found in pQ, specified by = (—1,0,0,0) 
and 

To [i ^ Ho . , 3 p 7 5 2 
^-^^l + ^smh p,F,^-,-,-,-smh 



(62) 



When we wish to investigate numerical properties of this solution, or of generalizations of 
it, we will usually focus on the parameter choices 



To = 5.55 



Ho = 0.33 , 



(63) 



These parameters, along with 



1 



4.3 fm 



(64) 



were identified in [T] as providing an approximate match to head-on gold-gold collisions at 
top RHIC energies, a/satw" = 200 GeV. 
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4.4 Properties of the viscous solutions 



Note that based on (60), since Fd is symmetric in x, T(p) is negative for sufficiently early de 
Sitter times. U Tq/Hq <^ 1 we can approximate Fd by a power law for very negative values 
of its argument and find the critical time, p*, for which T{p^) = 0. We find 



d-l 




. (65) 



For (i = 4, and for the parameter choices indicated in (63), which according to the analysis 
of [1] are similar to the values of the QGP generated at RHIC, we find that according to (65) 
p* ~ —6.09 which is fairly close to the numerical solution of T(p*) = which is p* ~ —6.15. 
At late de Sitter times the de Sitter temperature approaches a constant value of 

f (66) 

Negative temperatures is a pathology which signals that the hydrodynamic approximation 
should not be used at these early times. In physical terms, the shear is very strong at early 
times, so the viscous correction terms in the Navier-Stokes equations become comparable 
to or larger than the pressure terms. Failure of the hydrodynamic approximation can be 
characterized by the Knudsen number Kn which is given by the ratio of the mean free path 
to the typical representative scale of the flow. If the Knudsen number is much larger than 
one it implies that the continuum assumption necessary for the hydrodynamic approximation 
has broken down. In a conformal theory the mean free path is proportional to the inverse 
temperature. For the SO{d)-fiow in de Sitter space we use 



Kn=i 
T 



r 
f 



(67) 

As expected, the Knudsen number diverges at p = p*. It vanishes at p = 0, in the p — )■ oo 



limit, and at the zero of T'. A careful analysis of (60) shows us that such a zero exists only 
when p is positive and when the viscosity is large. 



Ho'^^W '''' 
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Figure 1: A plot of the Knudsen number defined in (67) for the viscous 5*0(3) symmetric 



flow given in (60). The Thick red line marks the location of the expected instability at early 



times discussed in the text. 



For = 4 we find that (68) reduces to 

To , v^r(-l/3) 



Ho 



< 



2 r(i/6) 



0.65, 



(69) 



A plot of the Knudsen number for the physically interesting values of d = 4, Tq = 5.55 and 
Hq = 0.33 is given in figure [1} 



4.5 Second order hydrodynamics 



In (28) we have written down the most general hydrodynamic stress-energy tensor, consistent 
with conformal symmetries up to first order in gradients. As is well known first order 
relativistic viscous hydrodynamics suffers from a causality problem (see, for example, [15]). 
That is to say, at large enough momentum superluminal modes may be present in the 
hydrodynamic solution. This isn't a problem per se since solutions with large momenta are 
outside the regime of validity of the hydrodynamic approximation. However, these modes do 



pose a problem when constructing a numerical solution to the Navier Stokes equation (32). 

In [To] the causality problem has been addressed by considering contributions to the 
stress tensor which contain two derivative terms. It can be shown that in the absence of a 
charged current and in a conformally flat spacetime the energy momentum tensor takes the 
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form 



d-1 



(70) 



1=0 



where . . . denote expressions involving three derivatives of the hydrodynamic variables, Aj 
are second order transport coefficients which scale like the temperature squared and the S*^*) 
are given by 



S(o) 
2(1) 



A ^(2) _ , ,A 



0'(mA^ u) 



S(3) 



(71) 



where 



- -P^P: {d,u, - d^u,) . (72) 

See [ini [12] for details^ 

We would now like to find an 5*0(3) x 5*0(1, 1) symmetric solution to the energy con- 



servation equation (32) in the presence of the second order hydrodynamic corrections given 



in (70). Assuming that = and restricting ourselves to ci = 4, we find that the energy 



conservation equation reads 

4 2 4 2 

a' H — a tanh p HoA/a(tanh p)^ H — Li tanh p H — L2 tanh p (sech p)^ = (73) 

3 3 3 3 



where y/a^ = and 



e'/'L, = 3 (Ao - Ai) 



ei/% = -(Ao + 2Ai) 



(74) 



The transport coefficients A2 and A3 do not enter into ( 73 ) since the vorticity Uajs vanishes 



for the velocity field in (47). Unfortunately, (73) is non-linear and difficult to solve. In the 



case were Li = L2 = it reduces to (55) with (58) whose solution is (60). On the other 



hand, if we set Hq = but keep Li and L2 non zero we find that 



a 



rp2 

-'o 



(cOshp)4/3 



(cosh pY 



(75) 



is a solution to (73). One can check that when converting this fiow back to R ' and taking 



the g — )■ 00 limit with T^q fixed, then the standard solution to Bjorken fiow with second 
order viscous corrections HOl is recovered. 



^When comparing to [10] it is useful to note that our Aq is related to Ttt through r/r^ — Aq. 
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p 



Figure 2: A plot of the Knudsen number defined in (67) for the 5*0(3) symmetric flow with 
second order corrections given in (75) with Li and L2 as in ( [76] ) and Tq = 5.5. The Thick 
red hues mark the zeros of T. 



For A/" = 4 Super Yang-Mills theory with gauge group SU (N) we find that at large N 
and large t' Hooft coupling [121 HH] , 



N 



63/4 



TC 



N 



127ry6 



;i-ln2)>0 L<2 



N 



127r\/6 



(4-ln2)>0. (76) 



Since Li > one finds that the solution (75) has vanishing temperature at times p = ±p*, 
T{±p^:) = 0, where p^, is defined through 



Li(cosh2/=^ p,f- fo(cosh2/3 - L2 = . 



(77) 



The existence of a critical value of p for which the temperature vanishes is unusual at first 



sight, but we have discussed a similar situation in section 4.4 when we considered first order 
viscous corrections (Hq 7^ 0, Li = L2 = 0) to the 5*0(3) symmetric flow. When we tried 



to run the viscous solution (60) backwards in time, we reached a regime where the viscous 



terms were large relative to the non viscous ones and the hydrodynamic approximation broke 



down. A similar situation occurs for the second order solution ( 75 ) as can be seen from figure 



|2j What is strange about ( 75 ) is that due to the time reversal symmetry of ( 73 ) when setting 
Ho = there exist two critical values of p. This may be a peculiarity of fluid flow in de Sitter 
space, or it may indicate a breakdown of the hydrodynamic expansion. 
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4.6 Angular momentum in 2 + 1 dimensions 

The case d = 3 is special because the symmetry group 5*0(2) x S0{1, 1) x Z2 doesn't fully 
determine the velocity profile. As noted previously, it is consistent with this symmetry to 
have a non-zero velocity ug around the 5*^ part of ^5*2. In the (p, 6, rj) coordinate system, 
the general 5*0(2) x 50(1, 1) x Z2-invariant velocity profile is 

where 5 is a function only of p. The total momentum around 5^ per unit rapidity is 

Pg = - dO cosh p top , (79) 



and the 5*0(2) symmetry guarantees that Pg is constant as a function of de Sitter time p. 
In the inviscid case with zero chemical potential, one can easily show that 

A = 37r-^coshV, (80) 

and that besides dPg/dp = 0, the only other equation implied by conservation of T^j, is 

di) 1 — v"^ 



dp 2 — v' 



-{)tanhp. (81) 



The equation (81) is separable, and the general solution of physical interest is 

V = ^£(v/2T^-£), (82) 

where 

£ = 4sechp (83) 

and io is a constant of integration. When ^ 1, the fluid moves close to the speed of 
light over a large region of dS2 x R. We will always pick the positive branch of the square 



roots in (82), thereby obtaining solutions with v > 0; but by choosing the negative branch 



on the outer square root, one can find solutions with v < 0. Plugging (82) into (80), one 
immediately obtains 

„ Pg £3/2 1 _ e(^/2TP - e) 



37r 



e2 



5 Vv2 + £2_£ 



4) 
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Temperature 



I- 4 




rapidity 




Figure 3: The temperature and rapidity in the future wedge of R^'^ for the flow specified in 
dS2 X R by the formulas (82) and (84), with £o = 1 and units chosen so that q= 1- 



Taking the — > limit together with Pe/^-o 
momentum solution (IgTI). 



1/2 



— )■ constant one recovers the zero angular 



The formulas (46) (with d = ?>) can be used to recover the energy density and velocity 
profile in the future wedge of R^'^. The velocity field in the (r, x, r/) coordinates takes the 
form 

/I 1) \ 

:,0 , (85) 



v y I — J 

and the flow is completely specified by the two functions f (r, x) and e(r, x), whose dependence 
on r and x is entirely algebraic (involving rational functions and square roots only), but too 
complicated to reproduce here explicitly. In figure |3] we show an example with 1^ = 1^ 
plotting the temperature T = e^l"^ and the rapidity /3 = arctanhf. It is interesting to note 
that the local energy density is symmetric under x — )■ —x even though the rapidity is strongly 
asymmetric. 

It should be possible to construct flows similar to the one described here in dimensions 
where the sphere has vector fields which are everywhere non-zero. It may also be possible 
to incorporate effects of shear viscosity and/or a non-zero chemical potential. 



5 Embedding in AdS space 



In [12] it was shown that any solution to (32) with zero chemical potential corresponds to 



a black hole configuration in an asymptotically AdS^ geometry, provided the appropriate 
transport coefficients take on their holographic values. 

Explicitly, the construction of jl2] is as follows. Suppose that u^, and e solve the d = 
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4 Navier stokes equation to order 0{d'^). By (9(9^) we mean terms which include two 
derivatives of the velocity field and energy density. For example, V^e is 0{d'^), but M^Vo-e 
is 0{d). Then, the line element 



ds'^ = —2u^dx^{dr + rA^dx") + r'^ g ^^dx^ dx'' 



1 



u^Uudx'^dx'' + br^F{br)af,udx^dx'' (86) 



with 



and 



irT 



F{br) ~ ^ + 2 ^"^^^^^^i^"^) ~ 2 arctan(6r) + - In (l — (6r) ^) 



is a solution to the Einstein equations 



^7) 



(89) 



provided terms of order 0(5^) are neglected. An extension of the solution (86 ) which includes 
second order corrections and a non vanishing chemical potential can be found in [TBI [13 1131 
[18] . See [I9] for a review. 

We emphasize that any solution to the conformally invariant version of the relativistic 



Navier-Stokes equation (32) can be written as a black hole solution in an asymptotically 
AdS geometry (with a boundary theory on R^'^) by using (86). An explicit and well known 



example of a solution which has been rewritten as a black hole solution to (89) is Bjorken 
fiow [201 EH [221 123]. Going to Bjorken coordinates. 



ds = Qfj^iydx^dx^ = —dr + dxj_ + t dr] , 



(90) 



Bjorken fiow is the following solution to (32) 



u' 



Toy/' _Hq 

tJ 2t 



(91) 



where tq is an integration constant, rj = H^e/T (which is a slightly different relation than 



(58)), and all other components of vanish. The AdS black hole dual of Bjorken fiow can 
be obtained by inserting (91) into (86), with 1/6 = ttT and setting the shear viscosity to its 
holographic value, Hq = l/(37r). Unfortunately, it is not so straightforward to check that 
(86) with (91) is a solution to (89) to order O^d"^). The difficulty lies in keeping track of 
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covariant derivatives of b (proportional to the inverse temperature) and u'^. One way to keep 
track of these derivative terms is to rescale the r coordinate so that 



T = Qf. 



(92) 



Then (91) takes the form 



u 



Q 



T[b{T) 



1/3 



Hp 
2Qf 



with To an integration constant, and after inserting (93) into (86) and using 



(93) 



ds"^ = gf.^dx^dx" = -Q'^dr'^ + dx\ + Q'^r'^drf, 



(94) 



one can easily check that (89) are satisfied to order 0{Q ^) in a large Q expansion. 
For the S'0(3) symmetric solution, which we rewrite as 



up 



-Qt tT = (coshp)-2/3 (to + ^^{smhp)^F^ Q, 



3 7 5 

6' 2'' 



■(sinh p)' 



on R^'^ with line element 



(95) 



ds"^ = g^udx^dx" = Q^T^ {-dp^ + (coshp)^ [dO"^ + sin^ ed(i?)) + Q'^drf 



(96) 



where 



Tj = HqeT ^ e oc T'' r = (cosh p cosh 6i — sinh p) ^ , 



the procedure is similar: We insert (95) and (96) into (86) with h = l/(7rT) and Hp = l/37r 



(97) 



to obtain a solution to the Einstein equations (89) up to order 0{d'^) which is equivalent to 
0{Q~'^). The extra factors of Q relative to (40) provide a simple book-keeping device, much 



like (92), which we use to keep track of gradients of the thermodynamic variables. To derive 



the solution (95) with these extra factors of Q one can rescale the spacetime coordinates of 



the Minkowski metric by a factor of Q, similar to (92), and make appropriate changes in 



equations (37) to (38). 
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6 Anisotropies I: setup and equations of motion 



The preceding section makes clear that passing from the future wedge of R'^"^'^ to dSd^i x R 
through a Weyl rescahng is a useful trick for generating exact solutions of conformal relativis- 
tic hydrodynamics. In this section we want to explore its utility in treating perturbations 
around the 5*0(3) x 5*0(1, 1) x Z2-invariant flow with zero chemical potential. In principle we 
could include chemical potential and/or consider arbitrary dimensions. But the case d = 4 
with /i = is the most interesting for applications to heavy ion phenomenology, so we will 
focus entirely on this case. 

The calculations in this section are somewhat reminiscent of the analysis of anisotropies 
in expanding geometries, which form an important part of theoretical cosmology (see for ex- 
ample [23] for a pedagogical summary). There are two important differences. One is that we 
are interested ultimately in the contracting half of ^5*3 (more specifically, the Poincare patch 
parametrized by (r, 77, 0)), whereas cosmological perturbations are more often consid- 
ered in expanding geometries. The other is that gravitational fluctuations are not included; 
the curved geometry in which we work is nothing more than a fixed Weyl rescaling of the 
future wedge of Minkowski space. All we are doing is solving the equations of relativistic 
hydrodynamics on a fixed, curved background, with perturbations treated (mostly) at linear 
order. 

The analysis of hydrodynamic anisotropies in hydrodynamic flow of the quark-gluon 
plasma at RHIC has also been recently reconsidered as a mechanism for explaining two- 
particle correlation measurements: see for example [3], SI El Ej [71 [8] . We compare our classi- 



fication of hydrodynamic perturbations and those of [8] in section 8.3 



6.1 Hydrodynamic perturbations 

As usual we do most of our calculations in the dS^, x R conformal frame, with coordinates 
= {p,6,(f),ri), and we use a hat to denote quantities specific to that frame. Also, we 
use a subscript b to denote unperturbed background quantities. Thus, for example, the 
background temperature of the 50(3) x 50(1, 1) x Z2-invariant solution which we perturb 



around is denoted Tf, in this section, whereas in (62) it was denoted T 



Perturbed fields may be parametrized as follows: 

f = n{l + 6) = (-l,z/„z/,), (98) 

where 6, z/j, and z/^ are all to be treated at linear order, and i runs over the two 5*^ directions. 
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By Ui we mean the two S"^ components of the four- velocity. Note that hnear corrections to 
the temperature are Tf,6. Corrections to the Up component of velocity enter only at quadratic 
order, so we suppress them. 6, z/,, and can depend on all four coordinates (p, 9, (j), rj). Our 
main strategy for computing these perturbations will be to decompose them into scalars 
and vectors of the 50(3) isometry group of dS^. Explicitly, a basis for hydrodynamic 
perturbations is 

Mp,9,(P,r]) = u,ip)d,S{9,^)e'''-'^ + Mp)Vi{9,^)e'''-'^ (99) 
u,{p,9,<P,v) = '^viP)SiOA)e'''^\ 

where, by assumption, 

~g'^ViV, = 0, (100) 

and g^^ and Vj are the metric and the covariant derivative of the unit 5*^ parametrized by 
{9,(f)). Symmetry arguments guarantee that the equations of motion for z/^, are decoupled 
from the equations of motion for Ug, z/^ and 6. Thus, in solving for the vector modes Uy, we 
can set the scalar modes to zero and vice versa. 

Practitioners of hydrodynamics will not be surprised to learn that the dynamics of the 



vector modes is purely diffusive. In section |6.2[ we give a fairly complete account of the 
vector modes. The scalar modes (5, Ug, z/^) are more complicated, exhibiting a combination of 
diffusive and reactive behavior. The early time instability mentioned in the introduction is 
also found in the scalar sector. We explain the basic equations underlying the scalar modes 
in section 6^, and in section [7] we study some of the solutions to these equations. 



In order to get physical insight into what the anisotropies represent, we must transform 
from the dS^ x R conformal frame back to more standard coordinates on the future wedge 
of R'^'^. This is straightforward. The flat space form of the perturbed flows is easily found 
to be 



T 


= T,(l + 5) 




b 99 


Ut 






b 99 


u± 






= ru^ 


Ur, 





(101) 



where by T^, and we mean the unperturbed flow obtained in (60) and (47) transformed to 
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R'^'^ using (46). 



6.2 Equation of motion for the vector modes 

A basis of divergenceless vector-valued functions on 5*^ is provided by the vector spherical 
harmonics [25| 1261127] 

9im{9, <p) = og H — d0 , liU2j 

sm (J sm (J 

which satisfy 

V.<l>L = and W.<I>L = - [^(^ + 1) - 1] $L for £=1,2,..., (103) 

where indices i are raised and lowered with the metric cjij of the unit S"^. When £ = 1, $' is 
a Killing vector of S"^. 



Consider now the excitations (99) with all the scalar perturbations (5, Ug, z/^) set to zero. 



and with Vi = Qij^j^. The equation of motion for which follows from (32) takes the form 

z/;(p) = -T,{p)Mp) (104) 

where 

= - tanh^ p H — ^ ° , tanh^ p 

- — ^ [-64 + 18£{£ + 1) + 9A;J + (16 + 9/cJ) cosh 2p] sech^ p tanh p 

and primes denote derivatives with respect to p. 

6.3 Equation of motion for the scalar modes 

A basis of scalar functions on S"^ is provided by the spherical harmonics Yfm{9,(f)), which 
satisfy 



27 



Consider the excitations (99) with the vector perturbations z/^ set to zero, and with S = Yim- 
If we define 

/A 



w 



(107) 



then instead of a single equation of the form (104), energy-momentum conservation gives us 
the coupled equations 



w 



(10^ 



where F, is a 3 x 3 matrix with all entries non-zero: 



Hq tanh^ p 

l)sechV 



Hq tanh p — Th 



ikr, 



3Tb 

2Ho tanh p + T), 



3Tfe 

2Hq tanh p 

^ + 1 

Ho tanh p — 2Tb 

HoTfe [-4(3£(£ + 1) - 10) sech^ p - 9k'^ - 16] + GHg tanh^ p + ST^ tanh p 



6T. 



Hq tanh p — 2T), 



(109) 



rs,12 

rs,i3 
rs,2i 

rs,22 
rs,23 

rs,32 
rs,33 



Alert readers will notice that there are three possible ways for divergences in components of 
to arise. The first route to a divergence in components of Fg is for Tb to vanish. This 



2Ho tanh p - 4% 



3Ho tanh p + Tb 



Hq tanh p + Tb 
i£{£ + l)B.okr, sechV 



4 


Ho tanh p + Tb 








fb 


{9i{i + 1) - 4)Ho sech^ p + 4Ho [3A;2 + 4] - 


STb tanh p 


+ 12B.I tanh^ p 


12Tb 


Hq tanh p + Tb 





occurs only at the moment when Tf, first becomes positive. As discussed in 4.4 we don't 
trust hydrodynamics in this region. The second route to a divergence is for Hq tanh p — 2Tb 
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to vanish. Noting that 



TKp) = - tanh p (Hq tanh p — 2Tb) 



(110) 



we see that the vanishing of Hq tanh p — 2T is equivalent to a zero of away from p = 0. 



The third route to a divergence is for Hq tanh p + T;, to vanish. It is clear from ( 105 ) that T 



also diverges when Tf, = and when T^' = 0, but not when Hotanhp + = 0. In section [7| 
we will explore the physical consequences of these divergences. 



7 Anisotropics II: special limits and stability 



While linear, equations (108) and (104) cannot in general be solved analytically. One can 



either resort to numerics, or consider limiting values of the parameters which make the 
equations tractable. A key question is whether the 5*0(3) x 50(1, 1) x Z2-invariant flow is 
stable against small perturbations. The answer we will reach is that there are instabilities 
related to flow in the rapidity direction. In the viscous case, these instabilities arise in the 
linearized approximation only at early times — earlier than the times at which hydrodynamics 
is typically initialized in hydrodynamic simulations of heavy-ion collisions. 



We start in section 7.1 with the vector modes, and then continue in section 7.2 with the 



scalar modes, investigating special limits and stability. We provide in section [7i3 an explicit, 
analytical account of how these perturbations grow large and qualitatively change the flow. 



7.1 Vector modes 



Formally, the solution to (104) is 



ly^ip) = exp 



dpTy{p) \ v^{pi) 



'1111 



where pi is an arbitrary initial time. The integral in (111) cannot be done explicitly in 



the viscous case because r„(p) contains hypergeometric functions. But in the inviscid case, 
where 



one flnds immediately that 



= — tanh p , 
3 



Uy{p) = (coshp)2/^z/40) 



'1121 



(113) 
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The exponential growth of Uy at large positive and negative p seems to indicate departure 
from the regime of validity of linearized perturbation theory. But this is not the case. The 
correct validity criterion is that the velocity vector should remain uniformly non-relativistic. 



From (98) and (99) we see that the velocity vector is non-relativistic precisely if 



1 I ^ cosh p , 



(114) 



and this is evidently satisfied by (113). So there are no instabilities in the inviscid case, and 



since (113) has no dependence on i, the general linearized vector flow can be written down 



immediately as 



Uy{p,9,(t),r]) = {coshpf^^Uy{0,9,(f),r]). 



(115) 



Intuitively, the main effect of viscosity should be to cause perturbations to decay, relative 
to their behavior in the inviscid case, as one moves forward in de Sitter time. We will consider 
two limits of r„ that, for the most part, confirm this intuition. The first limit is a small 
viscosity, short wavelength limit: 

Hq small; £ large; fc^ large; f^Hg finite; /c^Hq finite; p finite . (116) 
The second one is a late time limit: 



Hq finite, non-zero; 



finite; 



kj^ finite; 



p > 1. 



(117) 



Evidently, these limits are non-overlapping, so they offer complementary insights into the 
effects of viscosity. 



In the small wavelength, low viscosity limit (116) 



2 3Ho(£2 + fc2coshV) 
— tanh p H ^ h 

3 4To (cosh p) 4/3 



:il8) 



The omitted terms in (118) are suppressed relative to the ones shown either by positive 



powers of Hq, or by negative powers of i or /c^, or both. Because the second term in (118) 



is everywhere positive, it indeed suppresses perturbations at late de Sitter times relative to 



the inviscid case. It so happens that the integral (111) can be done explicitly using the 
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approximation (118), with the result 

3£2 



i^v{p) = exp 



3Ho 
4fn 



(cosh py/'^ 

X (coshp)2/3z/^(0) 



1 1 3 
2' 6' 2' 



sinh p 



sinh p 



(119) 



When fc^ = 0, the exphcit expression (119) leads to the following results for the large p 
asymptotics of vector perturbations: 

m^Hq ^r(-i/3) 



e ^'0 



^(1/6) (coshp)2/V„(0) ^ e (nr^.h n^2/3 



(coshp)^/V^(0) for p > 1 



e r(i/6) (coshp)2/=^z/^(0) 



0.97^^ 



(120) 



e""' (coshp)^/V^(0) for p < -1 . 



It is evident from the form of (118) that the suppression of vector modes becomes much 



stronger when one makes fc^ non-zero. 



In the late time limit (117) 



1 e^p/^ + - + 

2^/'2fn + Ho^S(^ 3^ ^ 



:i21) 



r(i/6) 



From the first term one sees that the behavior of vector perturbations is entirely stable at 
late times provided 



Ho 



> 



^r(-i/3) 
2r(i/6) 



0.647. 



'1221 



Comparing (122) to (69), we see that (122) is precisely the criterion for to have no zeros 
at positive p. 



Now let's consider the case where the condition (122) is violated. Then for A;^ 7^ 0, Ty 



becomes exponentially negative at late de Sitter times. As a result, vector perturbations 
with ^ have rapid growth (proportional to the exponential of an exponential of p) at 



late de Sitter times. There is a further peculiarity in Ty when the condition (122) fails: Ty 



has a pole at a finite positive value of p, call it pi, where T' vanishes. Let's consider the 
effects of such a pole in some general terms, starting by assuming the form 



r.(p) 



7 

p- pi 



finite . 



(123) 
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From (111) one immediately concludes that 



Uy{p) ^ Ky\p - pi\ 



(124) 



for p close to pi . Thus if 7 > there is a divergence in , signaling the failure of the linearized 
approximation; while if 7 < there is no divergence, only some failure of differentiability 
which is presumably smoothed out in the full non-linear theory. Inspection of shows that 



the pole that arises when (122) fails has 7 < 0. Thus there is no problem with the linearized 



approximation until the strongly negative behavior of the first term of (121) drives vector 



perturbations rapidly out of the linear regime. We note that for the physically interesting 



parameters (63), the left hand side of (122) is about 26 times larger than the right hand 



side. So fluid flows relevant for heavy-ion phenomenology seem to be safely away from the 
instabilities that we have described in this paragraph. If viscous corrections are so large that 



(122) is violated, one might reasonably ask whether fluid dynamics is applicable at all. 



7.2 Scalar modes 

The scalar sector is considerably more complicated due to the mixing among S, Ug, and z/^. 



The plan of this section is to investigate solutions to (108) and stability of the scalar sector 



as far as possible in parallel to our study of vector perturbations in section 7.1 In particular, 
we will start with the inviscid case (with one additional simplification), then explain viscous 
damping in a small wavelength limit, and finally examine stability by studying the late time 
limit and behavior near singularities of components of Tg. 

The inviscid case is not analytically tractable in general due to mixing among 6, Ug, and 
Uri- This mixing goes away when fc^ = 0, in which case (with Hq also set to 0) we find 



1 



h 1) sech^ p 
tanh p 







tanh 



(125) 



Evidently, decouples and obeys the same equation as the vector modes. So 



z/^(p) = (coshp)2/3z/^(0). 



(126) 



In the case of the vector modes, we remarked that the exponential growth in i^vip) indicated 



in (113) at late times is not a problem, because the condition for flow to be non-relativistic 
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is {uyl ^ cosh p. The situation is less favorable for i/„: the condition for the flow to be 



non-relativistic is |z/„| ^ 1, which evidently is violated by the solution (126) for any non-zero 



^'r,(0). When z/^(p) gets big, the linear approximation is not justified. We will soon see that 
viscosity cures this particular instabihty, so it can be regarded as an artifact of setting the 
shear viscosity strictly to zero. 



It is straightforward to show starting from (125) that Us can be eliminated through the 
equation 



+ 1 



cosh p 



dS 
dp 



(127) 



and that S obeys the simple second order equation 



d"^ 4 ^ d 



sech^ p 



0. 



(128) 



This last equation can be solved in terms of associated Legendre functions: 



5{p) = (sechp) 



2/3 



K^P 



2/3 , 



2/3 



tanh p) + K+P '[ 



+ |^l+12£(f+l) 



(tanh p) 



(129) 



where K± are integration constants. If only is non-zero, then Vs remains finite as p — t- — oo 
but grows as e^^/^ as p — ?■ +oo. As with vector modes, this late-time growth is not a concern 
because the criterion for departure from the non-relativistic regime where linear perturbation 
theory is good is \vs\ <^ cosh p. 

A good way to understand the overall behavior of the scalar modes is to consider the 
eigenvalues of Tg- When Hq = A;^ = 0, these eigenvalues may be read off easily from (125): 



1 



1 



A-t = — tanh p ± - sech p \ sinh p — 3£(£ + 1 



A. 



(130) 



tanh p . 



Clearly, A^ is the eigenvalue related to the diffusive mode f^, while \± relate to the (5, Vs) 
modes. These modes are oscillatory when \± are complex, which is to say when |p| is not 
too large. At large |p|, the {S,ris) modes are non-oscillatory. 



Now let's consider what happens when shear viscosity is added. As in section 7.1 



find it useful to consider two complementary limits. The first one is a small viscosity, short 
wavelength limit: 



Ho small; 



large; 



0; 



^Hq finite; p finite . 



[131] 
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The only reason to require A;^ = as part of the hmit (131 ) is that it substantially simplifies 
the equations. The second limit we consider is 



Hn finite, non-zero; 



finite; fc^ finite; p ^ 1 , 



(132) 



which is identical to (117). 



In the small wavelength, low viscosity limit (131) 



1 

+ 



Hn 



To(coshp)4/3 



+ 1) sech^ p 
^ tanh p 


/o 









I tanh p j 
?+ l)tanhp 

e 





3/)2 



+ 



(133) 



where, as in (118), the omitted terms are suppressed relative to the ones shown either by 



positive powers of Hq, or by negative powers of ^, or both. Clearly, the matrix in the second 

3^^Ho 



line of (133) shifts A.^ by + 



4To(coshp)4/3 ■ 



just as in (118). So once again, z/^ formally behaves 



just like the vector modes. The (5, modes are more complicated. Numerical investigation 



shows that the shift of \± induced by viscous corrections in the limit (132) is not uniformly 



positive when \± are purely real; but it is usually positive. In particular, one can show by 
explicit calculation that during the oscillatory phase, where \± are complex, the effect of 
viscous corrections on Re \± is to shift it by + 



2f ( cosh°p)4/» • ^ ^^^^ generally in line with 
the intuitive expectation that the viscosity damps out all perturbations, relative to their 
inviscid evolution, as p becomes more positive. 



Now let us turn to the late time limit (132), where one finds 





3 '^r? 





\ I 

" 3 



+ 



1 ^ 

3 



-k \ 

3 n 





+ 0(e-2p/3) 



(134) 



where 



Hi 



Hr 



25/3 9^ , TT v^r(-l/3) • 

2Jo + ilo p(i/g) 



(135) 



Note that there is no £ dependence in V ^ to the order shown in (134). The eigenvalues of F 
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as given by (134) are easy to find explicitly: 



Ai = ^(4 + fcJ-^4-31fc2 + A;4^ 
A2 = ^ (4 + fcj + ^4-31fc2 + fc4) (136) 

A3 = niky^/' + 1 . 

It is clear from the first term in A3 that late stability for k„ ^ holds only Hi > 0, which 



is the same condition as the bound (122). The pathologies that arise when this bound is 
violated, together with the observation that violation of the bound is far from the regime 
of interesting parameters in heavy ion collisions, dissuade us from considering Hi < here. 
Instead we will focus on the case Hi > 0: then all three eigenvalues Aj have positive real 
parts for arbitrary real krj. This is generally indicative of stability. Numerical exploration 
confirms that all three quantities {6, Us, z/^) are driven to small values at late times p. Thus, 
shear viscosity cures the late de Sitter time instability in z/^. 

The only remaining route to instability is for some component of Fg to diverge at finite p, 
in the region where Tb > 0. This indeed happens, in two different ways. The first way is that 
if Hi < 0, three components of Fs diverge at some pi > 0. As stated earlier, this instability 
is associated with zeros of T^', and we will not explore it further. The last type of divergence 
arises when Hgtanh p + TJ, vanishes. This occurs at some de Sitter time p2 which is negative, 
but less negative than the time p^, where Tf, first becomes positive. Using the identity (110), 



one can show that the Knudsen number at p = p2, as defined in (67), is I/Hq. Keeping in 
mind that Hq = 0.33 is in the ballpark of realistic numbers for a heavy ion collision at top 
RHIC energies, we see that fiuid dynamics is not particularly reliable at p = p2. Indeed, 



p = P2 corresponds to proper time r ~ 0.05 fm/c at x_l = (see section 8.2 for a fuller 
discussion) which is a substantially earlier time than is usually considered reasonable for 
hydrodynamic simulations to start. Despite these observations, we wish to understand more 
fully what happens to the hydrodynamic description at p = p2. 

To simplify our investigation of the divergence at p = p2, let's start by assuming /c^ = 0. 
Then i/^ decouples from 6 and Us, and the only component of Fg that diverges at p = p2 is 
rs 33. Thus we only need to consider the evolution of z/^. Near p = p2, one can show from 



(109) and (110) that 



T 3 
33 = — — + finite with 7 = 1 + -i(i + 1) sech^ ps . (137) 
p - P2 4 
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Evidently, 7 > for all £. This means that z/^ diverges as \p — P2\~'^ ■ The divergence is 
stronger for larger I. Numerical exploration indicates that when A;^ 7^ 0, the divergence 
is still present; in fact it spreads to all three quantities (5, z/g,!/^). To say more about this 
divergence, we need to go beyond the linearized approximation. We take some steps in this 
direction in section 17.31 



7.3 Non-linear treatment of rapidity perturbations 



We learned in section (7.2) that the perturbation that most often exhibits instabilities is 
Vri- We would like to go beyond the linearized treatment of this perturbation. To do so, we 
assume £ = and fc^ = 0: that is, 50(3) x 5*0(1, 1) symmetry is preserved. But we allow 
the Z2 symmetry that sends 7] — )■ —r] to be broken. The most general hydrodynamic ansatz 
consistent with these symmetries is 



u 



p 



coshx(p), ue = u^ = 0, M^ = sinhx(p). (138) 



Temperature T (or, equivalently, energy density e) must also be a function only of p. 



With the ansatz (138), the Euler equations read 

dx sinh 2x 



dp 2 + cosh 2x 
d log T 2 cosh^ X 



tanh p 
tanh p 



(139) 



dp 2 + cosh 2x 

These two equations are readily solved because the first of them is separable. The general 
solution can be given in implicit form as 

cosh^ p = Ki sinh^ y tanh v 

(140) 

T = K2 csch X 

where Ki and K2 are constants of integration, which must have the same sign in order for 
T to be positive. We assume that Ki and K2 are positive, which is the same as choosing x 
to be positive. Note that at late times, 

Y ~ P log Ki 

^ 2 ^ ' (141) 

f ^ 2J^iK2e~''. 
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This is a non-linear version, for £ = 0, of the instabihty of the strict inviscid hmit which we 



noted following (126). 



Now let us pass to the viscous case. Plugging the ansatz (138) into the viscous Navier- 



Stokes equations (with vanishing chemical potential) leads to the following two equations: 



dPr, 

dp 



Pn 



— f cosh p T 

dp V y 167rHoT2 



sech p csch x sech x + 



3P2 



2567r2HnT6 



sech p csch x sech x ■ 



(142) 



where 



Pr. 



An cosh^ p T\ 
167r 



cosh p cosh X sinh x T 



1 + 



Ho 
f 



tanh p cosh x 



— smh X 
dp 



(143) 



is the momentum in the rj direction per unit rapidity. We were unable to find the general 
solution to the equations (142). However, an interesting special case is = 0, corresponding 



to flows that have no net momentum in the 1] direction. One way to achieve = is to set 
x{p) = everywhere. Then one is led back to the 5*0(3) x 5*0(1, 1) x Z2-invariant solution 
(62). Let us instead assume x 7^ 0, but nevertheless P^ = 0. Then the second equation 



in (142) can be used to conclude that T oc sech p. Also, the quantity in square brackets in 



(143) must vanish. Combining these constraints, one sees that the general solution to (142) 



with X 7^ and P^ = is 



cosh X = Ks cosh(p — po) 



-Ho-ft's sinh po sech p , 



(144) 



where and po are constants of integration. If = 1, and we assume that x is increasing 
as a function of p, then we have the simpler form 



x{p) = p-po T(p) = -Ho sinh Po sech p. 



(145) 



It is interesting to note some properties of the flows (144) in other frames. In fiat space, 

T 



-2gHo-ft'3sinhpo 



v± 



V^l + 2g2(r2 + xl) + q^r^ - x^^f 
u± 2q'^TX± 



(146) 



1 + q^T^ + q^x 



'2 rf'^ 
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The expression for v± is unchanged from the S'0(3) x S'0(1, 1) x Z2-invariant flow: it is 
simply the expression of the condition = in flat space coordinates. An expression for 
t>^ = can also be given, but it is long and unenlightening. It is worth noting, however. 



that for the special solutions (145), the range of is the whole interval (—1, 1) allowed by 
causality. 



In the AdS2 x S*^ conformal frame described in (41), the solutions (144) satisfy 



T = cosh pT = Tq 



Hn 



(147) 



where Tq = —HqK^ sinh po, and we use rings, as on T, to denote quantities in the AdS2 x S'^ 



frame. The simplicity of the formulas (147) suggests that further study of solutions with 



non-zero x would be easiest in the AdS2 x S"^ frame. 



Now let's compare the non-linear flows (145) to the linearized analysis of rapidity fluctua- 



tions. Numerical exploration indicates that if one starts at p < p2 with a slight perturbation 



of the 5*0(3) X 5*0(1,1) x Z2-symmetric flow (62) toward non-zero f^, with 



0, 



the resulting flow quickly converges approximately to (145) with po ~ P2- See for example 



figure |4j This convergence is not surprising: the initial conditions clearly specify a flow with 
small P^, and we know from the linearized analysis that x cannot stay small. So the general 
solution (144) is a natural candidate to describe the evolution. Picking po = p2 is natural 



because this is the time when the rapidity starts departing from the linear regime. It is less 
clear to us why the choice Ki = 1 works as well as it does in describing the flow for p > p2. 
Three additional points are worth noting: 



The flows (|144|) do not in general have a small Knudsen number. Using the definition 



and, for the special case 



(67), we have immediately Kn = which is small in some neighborhood of p 

\X'{P)\ 

f{p) 



and large elsewhere. It is perhaps fairer to define Kn 
Ki = 1, this leads to the expression Kn = which has similar qualitative features 

to the previous expression provided Tq ^ 1. 



Solutions to the non-linear equations preserving 50(3) but not 50(1, 1) can prob- 
ably be generated starting from perturbations with non-zero krj. Such solutions are 
significantly harder to study because the differential equations to be solved have two 
independent variables, r and t]. Presumably the non-linear evolution leads to rapidity 
gaps. Such gaps are not observed in heavy-ion data as far as we are aware. 



According to (137), there are modes with £ 7^ which show similar instabilities to the 
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X 



;r(-5.2) = 0.1, 
linear 



;r(-5.2) = o.i 



^-(-5.2)= 0.1, 
non-linear 




Figure 4: Left: The rapidity x as a function of de Sitter time p, both in the hnear ap- 
proximation (green) and in the full non-linear theory (red). Right: The temperature of 
the S'0(3) X S'0(1, 1) x Z2-invariant flow (green) and of a flow with non-zero x (red). The 
viscosity parameter is Hq = 0.33, as in (63). The specific solutions plotted are based on 
choosing initial conditions at pi = —5.2 such that T matches the flow (62) with parame- 
ters chosen as in (63), and such that x(~5.2) = 0.1, with dx/dp chosen in the non-linear 
treatment to match its value in the linearized treatment. The dashed black curves show the 



analytical solution (145) with po chosen so that T(0) matches exactly between the numerical 



and analytical solutions. 



one in aX = £ = which we have analyzed in this section at the non-linear level. 
It therefore seems likely that the full 3 + 1-dimensional dynamics is chaotic, with many 
sound mode instabilities, all related to the one studied here, but with non-zero i and 
fc^, developing simultaneously. 



Let us end this section with a discussion of singularities. The singularity at p 



P2 



noted in the linear theory at the end of section |7.2| is smoothed out in the non-linear theory. 
However, the non-linear theory has singularities of a new sort: If < 1 in ( 144[ ), then 
X — )■ at two finite values of p, call them p±, with p_ < p+. For p G {p^,p^), x{p) is 
undefined, and as one approaches p_ from below, x ~ VP- ~ P- Thus the flow becomes 
non-differentiable at p = p^. There is a similar singularity at p = p+. 



8 Applications to heavy-ion physics 



As already noted in [Ij and around equations (63) and ( [64j ), the S0{3) x 5*0(1,1) x Z2- 
symmetric flow with Tq = 5.55, Hq = 0.33, and q = 1/(4.3 fm) is in the correct ballpark to 
describe the hydrodynamic evolution of head-on gold-gold collisions at top RHIC energies, 
y/sj^j^ = 200 GeV. This description is obviously imperfect: neither the initial conditions 
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nor the hydro dynamic behavior are truly conformally symmetric in real heavy-ion collisions. 
Nevertheless, the 5*0(3) x 5*0(1, 1) x Z2-symmetric flow is useful because it correctly captures 
aspects of real heavy-ion collisions in a formalism which is to a large extent analytically 
tractable. The purpose of this section is to provide some translations between aspects of 
heavy-ion dynamics and the mathematical formalism based on conformal symmetry. In 



particular, we will map out in section 8.1 the approximate locations of thermalization and 



hadronization in the dS^ x R conformal frame; we will provide in section 8.2 several examples 



of perturbed solutions exhibiting elliptic and triangular flow; and we will explain in section[8T3 



how the natural treatment of initial conditions based on conformal symmetry compares to 
the cumulant expansion of [8]. 

8.1 The history of a heavy ion coUision in the de Sitter conformal 
frame 

In this section we aim to show where in dS^ x R thermalization and hadronization occur. We 
will continue to focus on head-on collisions at top RHIC energies (that is, y^s^N = 200 GeV 
for gold-gold collisions). Our results are presented in figure |5| which shows how the anatomy 
of the collision in the {t,x±) plane gets repackaged in the coordinates {p,0). 

First let's explain how we plotted the region where thermalization occurs in the {p,6) 
coordinates. Typical choices of thermalization time in hydrodynamic simulations include 
Bjorken r between 0.6fm/c and Ifm/c, though values as small as 0.3fm/c are probably 



acceptable. From (11) one can derive the formulas 



sechp sin^ , . . 

Using these formulas one can find the locus of points (p, 6) that corresponds to a given r (or 
to a given x±). In figure [s] we show the curves r = 0.3fm/c, 0.6fm/c, and Ifm/c. 

Hadronization occurs not on a specific time-slice, but through processes at characteristic 
temperatures. Chemical freeze-out, determined by thermal fits to hadron yields, occurs at 
T'chem ~ 170 MeV (see for example |2S]), whereas kinetic freeze-out, determined by thermal 
fits to the momentum dependence of particle yields of a fixed species, occurs at Tkin ~ 
110 MeV (see for example [29]). These temperatures are low enough that deviations from 
the conformal equation of state, e = (^^T^ are large, and conformal methods are probably 
no longer useful. Nevertheless, in order to get an approximate idea of where in dS^, x R the 



40 




Figure 5: (Color online.) An approximate history of a head-on heavy ion collision as seen 
in the dS^ x R conformal frame. Fluid dynamics can be used in the orange region labeled 
"QGP." If thermalization occurs earlier than r = Ifm/c, then fluid dynamics is also valid 
across part or all of the green region labeled "Thermalization." The thick curve labeled 
"total overlap" shows a characteristic time r = 0.07 fm/c for total overlap of the colliding 
nuclei. This curve ends at x_l = 7fm because that is approximately the full transverse radius 
of a gold nucleus. The line labeled "z/^ instability" shows where the instability described in 
sections 7.2 and 7.3 occurs. The parameters (63) and (64) were used in drawing all curves 



in the figure. We also used the S0{3) x 5*0(1,1) x Z2-invariant solution (62) in order to 
calculate the position of the curves depicting freezeout, and also of the lines showing where 
T vanishes and where the instability occurs. 
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validity of fluid dynamics ends, we will use the relation 




(149) 



with = 11, in order to find the locus of points (p, 6) where a given flow attains a definite 
Minkowski space temperature. The value f^, ~ 11 can be extracted from lattice data |30| |3T] 
for 1.2T, < T < 2Tc. In figure [5] we show the curves T = 110 MeV and T = 170 MeV. 



8.2 Examples of perturbed flows 



In sections ( 7.1[ ) and (7.2) we constructed analytic solutions for the linearized perturbations 
around the inviscid S'0(3)-invariant flow. When considering the S'0(3)-invariant flow as a 
toy model for the thermodynamic phase of heavy ion collisions, the £ = 2 perturbations may 
be thought of as elliptic flow and the i = 3 perturbations as triangular flow. 



Consider first the scalar modes (129) in the inviscid theory, with fc^ set to (i.e. the 
perturbation doesn't break S'0(1,1)). Let us also set = 0, so that the temperature 
perturbation is 

*(rt = A-_(sechrt=/'P^^'^^^j-^(-ta„hp)^ (150) 

The point of setting Kj^ = is that then the velocity perturbation Us is as small as it can 
be in the early time limit. An intriguing feature of these modes is that one can predict their 
amplitude at large p through a very simple formula: 
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lim 6{p) ~ ^/3 
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' = 2 


1.15 


for I 


' = 3 
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for I 
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-1 


for I 


' = b 



(151) 



The low odd multipoles are thus significantly enhanced relative to the low even multipoles, 
as measured by their propagation from p = — ootop = oo. It might seem from this result 
that we are predicting an enhancement of odd multipoles over even in heavy ion collisions. 
Unfortunately, this is not true, for two reasons. First, as is clear from figure [5} freezeout 
surfaces are far from being located at large p: indeed, with the assumptions that went into 
this figure, the kinetic freezeout surface only barely extends to p > 0. And second, viscous 



corrections are significant for our usual, quasi-realistic choice of parameters (63). The effect 
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viscous 1=2 



inviscid 1=2 




Figure 6: (Color online) Inviscid (solid) and viscous (dashed) perturbations of the 5*0 (3)- 



invariant flow. The temperature perturbations S in the inviscid case are given by (150) 



and the velocity perturbations Ug follow from (127). The viscous perturbations are obtained 
numerically, with initial conditions chosen at p = —1.97 (corresponding at = to r = 
0.6fm/c) to agree with the inviscid perturbations. 



of viscous corrections can be seen by numerically integrating the viscous equations (108), still 
with k„ = 0, and with initial conditions chosen at some reasonable initial time to be as close 



as possible to the inviscid form (151). In figure 6 we show the results of such integrations. 



The effort involved in finding the viscous curves in figure |6] is slight compared with full 
viscous hydrodynamics codes. 

In order to gain further intuition about the perturbed conformal flows, let us now study 
them directly in the future wedge of R^'^. We will not address hadronization, and for 
simplicity we will work in the inviscid theory; thus our results are not wholly realistic, and 
are meant mostly to convey in familiar coordinates what the perturbed conformal flows 
mean. We will focus on the low multipoles of the scalar modes, namely 



6''^ = (sechp)2/3p2/3 /-(-tanhp) 
6^'^ = (sechp)2/3p2/3 _(-tanhp) 
6^'^ = (sechp)2/3p2/3^^^^(-tanhp) 



5=^'^ = (sechp)2/3p2/3 ^ (_tanhp) 
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3,31 



(152) 

with corresponding perturbations of the velocity field (127). Our conventions for spherical 
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harmonics are given in appendix [bJ As mentioned earlier, because 5^'^ involves modes with 
m = ±2, it should relate to elliptic flow. We chose the particular admixture of ¥2^(6, (p) so 
that 5^'^ as a whole would be invariant under an 5*0(2) subgroup of 5*0(3) which contains 
rotations of 5^ orthogonal to the one generated by c?^. Both 5^'^ and 5^'^ are i = 2 modes, so 
their p dependence is the same. However, since 6 as well as p depends on r, the r dependence 
of 6^'^ and 5^'^ are rather different. The same reasoning applies to the 6^'^ and 6^'^ modes. 

Using ( |TT| ) one can readily express the perturbations (152) as functions of r, x\ and 
x^. In flgures [t] and [s] we show the energy density of the background flow ^ perturbed 
by the perturbations in ( |152 ). More precisely, we set T = (l + for i and m 

as in (152), and so the energy density plotted is (^1 + K^6^'"^Y with as in (18) and 
= 3/5, q = (4.3fm)~^ and eo = 880. We have used = 3/5 in order to make the 
main features of the perturbation obvious. The value for e has been taken from [Ij. While 
the temperature remains positive, the perturbations as plotted can hardly be described as 
small, so the linearized approximation may not be justifled for such a large value of K^. For 
tracking the time evolution of the energy density, we found it useful to plot the quantity 

'^^'^ ^ - e,{r,r,0,0) " e,(r, r, 0, 0) + " ^'^^^ 

Put in words, e is the perturbed energy density divided by the unperturbed energy density 
evaluated at = r with = 0. This quantity is convenient because its global maximum 
on any slice of constant Bjorken time r is close to 1. 

For illustrative purposes, we've also plotted, in flgure|9| vector perturbations of the form 
Ui = KvV^'^ with 

V''^ = (coshp)2/3-L ($3 _3 _ $3 3) (154) 

and = 3/5. As expected, the vector perturbations don't affect the rest frame energy 
density e and exhibit a non trivial velocity fleld. In the lab frame, the energy density is time 
dependent as is the energy flux. 



8.3 Parameterizing anisotropies 

Recently, the authors of [3j have provided an explanation of certain features of two-particle 
correlations at RHIC which relies on the propagation of fluctuating initial conditions. Roughly, 
stochastic fluctuations in the initial phase of the collision evolve into anisotropic perturba- 
tions of the hydrodynamic flow which, after hadronization generate the peculiar "ridge" and 
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Figure 7: The energy and velocity field for the 5^'^ (top) and b^'^ (bottom) perturbations 



described in equation (152). The energy density is lowest in blue regions and is close to its 



maximum in red regions. The two left columns show density plots of the energy density 
overlaid with a vector plot of the velocity field whose ri component vanishes. The rightmost 



plots are contour plots of e described in ( 153 ) and provide an illustration of the time evolution 
of the energy density. 
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Figure 8: The energy and velocity field for the 5^'^ (top) and 5^'^ (bottom) described in 



equation (152). The energy density is lowest in blue regions and is close to its maximum 
in red regions. The two left columns show density plots of the energy density overlaid with 
a vector plot of the velocity field whose ?7 component vanishes. The rightmost plots are 



contour plots of e described in (153) and provide an illustration of the time evolution of the 
energy density. 
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Figure 9: The energy and velocity field for the vector perturbation V^'"^ described in equation 



(154). The energy density is lowest in blue regions and is close to its maximum in red regions. 



"shoulder" observed in two particle correlations. In ^ this idea was formulated using a 
cumulant expansion of the early time perturbations of the entropy density which was then 
fed as an initial condition to linearized hydrodynamics and finally, using a Cooper-Frye 
hadronization scheme, to a final hadron distribution. 

It is interesting to compare our expansion of anisotropies in spherical harmonics with 
the cumulant expansion of [8jj^ To simplify the discussion, let's suppose that the stochastic 
anisotropies of the flow are due to anisotropies in the transverse energy density of the colliding 
ions which we approximate by a light-like shock wave: 

n = -V{x±)S{u) (155) 

where x± = and 

u = t-x^ v = t + x^. (156) 

The actual form of the transverse distribution f{x±) is model dependent. As explained in [T] 
(and anticipated in [HUES] on AdS/CFT grounds), the unique S'0(3)-invariant transverse 



This section is based on discussions with D. Teaney. 
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distribution is 



/50(3)(^±) 



(157) 



where E is the total energy of one shock wave. A more phenomenologically motivated choice 
is based on the Woods-Saxon profile: 



/ws(^±) 



E 



(158) 



-87ra3Li3(-e^/'^) J _^ ^ ^ ^yxiMx^-R]/a ' 
where E is again the total energy and Li„ denotes the polylogarithm function. The rms 



radius of the distributions (157) and (158) are the same provided 



Li3(-e^/") 
8a2Li5(-e^/'^) ' 



(159) 



Finally, following we can consider 



^Gaussian 



q^E 



TT 



2 2 



(160) 



whose rms width is the same as (157). It would be more precisely in the spirit of [8] for / to 



be the transverse density of wounded nucleons in a Glauber treatment of the early stages of 
a heavy ion collision. In central collisions, this density approximately follows the transverse 
energy density of one of the nucleons, so the distinction between our definition of / and the 
one used in [8] is not large. In what follows we will set g = -E = 1 for simplicity. 

In order to facilitate the comparison between the anisotropics parameterized by spherical 
harmonics in this work, and the cumulant expansion of [8j, it is useful to introduce complex 
coordinates and derivatives in the transverse plane (a;^,x^) 



x^ + ix^ 



d 



d 1 d 
+ 



dx^ i dx^ 



d 



X 

1 
2 



IX 



d 1 d 



dx^ i dx"^ 



(161) 



If one takes de Sitter time p — t- — oo, corresponding to Bjorken time r — > from above, then 



starting from the second equation in (11) one can show that the transverse plane is mapped 
to S"^ as follows: 



z = e'^iim- z = e-'^idJi-. (162) 
2 2 

This is the standard stereographic map, which itself is conformal. To find the Weyl factor Vt 



e 
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associated with the stereographic map (162), we note that the metric on S"^ takes the form 

idzdz 



d.r = de^ + sin^ 9 , 



(163) 



from which we conclude that 



Q(z, z) 



l + zz 



(164) 



Consider the 50 (3) -invariant transverse distribution (157) in complex coordinates 



SO{?.){Z,Z) 



1-3 



47r 



(165) 



We see that / maps to a constant on 5^ under the stereographic projection (162) if it 



is assumed to have conformal weight 3. So for functions /(z, z) not too different from 
fso{?,){.z,z), it is sensible to expand 



where 



oo I 

f{z,z) = n{z,z)-^Yl Yl /^.-^^ 

e=o m=-l 



,m\Z, Z) , 



(166) 



(167) 



with the coordinate mappings as defined in ( 162[ ). It is straightforward to show that for 
m > 0, 

(168) 



Xp,miz, Z) 



;Pe..m\zz 



(1 + zzY 

where Pe^miO is a polynomial of degree i — m with all its zeros at positive values of ^. 
Because Xg^miz, z) = ±X£^_m{z, z)*, there is essentially no new information in the Xi^^ with 
m < 0. We tabulate the first few Xi^m with m > in Table |2] 

The coefficients fe^m can be interpreted as moments of slightly unusual quantities. Defin- 
ing 

, ^^ f d'^X± f(x±)A(x±) 

for an arbitrary function A{x±) with good smoothness and fall-off properties, we see thalj^ 



{n{z,z)Xi^rn{z,Zy) 



f d^x± /(x_l) 



(170) 



^It is straightforward to generalize (166) and (170) to cases where / transforms with an arbitrary weight 



a under the stereographic map: all that changes is the powers of 17. 
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Assuming that J (Px± f{x±) = 1 we find, with the conventions of appendix^ 




TT 



1,0 



IGvr 



(1 - zz) 



(171) 



These moments can be assumed to vanish, since non-zero fi^i corresponds to a displacement 
of the distribution in the transverse plane, and non-zero fi^ corresponds to an alteration 
of its overall size, which can be compensated for by changing q. Low-order moments which 
cannot be assumed to vanish include 



2,2 




Z^ 



TT \1 + ZZ 



/2,1 




15 

2-K \1 + zz 



3,3 




^3 



l + zz 



(172) 



where we have assumed /i^i = in order to simplify the expression for /2,i. 

Now let us give a slightly simplified treatment of the cumulant expansion of [8]. The 
main idea is that a general function (or at least an appropriately broad class of functions) 
can be represented as 

/(f^) = / ^^^k^■^^eWCk^) ^ (173) 



where W{k±) admits a Taylor series expansion. In practice, for the scope of work considered 
in [8], it is generally enough to consider 



kl 
4 



where SW{k±) is treated as small. Thus, to linear order 

1 d 



1 + SW 



i dx 



TT 



Another way of expressing (175) is 



(174) 



(175) 



IT 



n=0 m = -n 

m-\-n even 



where we define 



n — m — n-^m 



Qn,m{z,z) = e^'^'d 2 (9 2 e 



(176) 



:i77) 



Clearly, the Qm,n are polynomials in z and z. We tabulate the first few Qn,m with m > 
in Table [2| It is clear from the table that ^ is qualitatively similar to Q2£-\m\,m- The 
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up to a constant 


e = o 


e = 1 


£ = 2 


£ = 3 


m = 


1 


1— zz 


l-4zz+z2z2 


l-Qzz+gz^z^-z^z^ 


1+zz 


(1+ZZ)2 


{1+ZZ)3 


m = 1 




2 


z{l—zz) 


Z(l— + Z z j 




1 + ZZ 


(l+zz)2 


{1+Z2)3 


m = 2 








22(1-ZZ) 






(l+zz)2 


{1+ZZ)3 


m = 3 








z3 
(1+Z2)3 



up to a constant 


n = n = 1 n = 2 n = 


3 


n = A 


n = 5 


171 = 


1 I- zz 




2 - Azz + z^z^ 




m = 1 


z z{2 - 


zz) 




z{6 — 6zz + z^z"^) 


m = 2 






z^{3 - zz) 




m = 3 


z^ 






z^{A - zz) 



Table 2: The first few ^ and Qn,m witli overall constant prefactors omitted. 

main difference is the powers of 1 + zz in the denominators of the Xi^m, which renders them 
uniformly bounded across the plane. The unboundedness of the Qn,m was a problem in [8], 
solved by an ad hoc regularization. 



The constants Wn,m in (175) are closely related to the moments Wn,m, W^^, and 
of [8] Fl To work out expressions for Wn,m in terms of moments, it helps to define 



kz = ^(,ki- ih) 



h = ^{h +ik2) 



(178) 



Then from (173)- (177) we read off 

oo n 

l + SW{k,,k,) = J2 E wr.,m{^k.)'^{tk,)'^ =e'^'^{e-''^^-''^') , (179) 

n=0 m^-n 

m + n even 

where we persist in treating 6W at linear order and we assume Wqq = 1. Matching terms in 



^One difference in notation between the present summary and [5] is that n and m have been switched: 
For us, m denotes the "magnetic quantum number" which indicates that a function contains a factor e*™"^, 
whereas in [8] n is the magnetic quantum number. 
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Taylor series expansions in kz and kz, one finds 



il-zz) 



(180) 



These moments are identical to those in (171 ), and we again assume that they vanish. Low- 
order non- vanishing moments include 



W3,3 



:i8r 



where we have assumed wi^i = in order to simplify the expression for w^^i. 



There are two reasons in heavy-ion phenomenology to regard the moments (172) as 



useful alternatives to the more commonly used ones in (181). First, they are based on 



5*0(3) group content, and the strength of the perturbations that they measure would be 
preserved by conformally invariant dynamics. The early time dynamics of a heavy ion 
collision is approximately conformal, and it seems natural to take advantage of this. Second, 
purely polynomial moments like (z^) have a strong bias toward surface effects. The moment 
{z^/ (1 + zz)'^) has less surface bias. 



9 Summary 

In this work we found an analytic solution for a viscous, conformal fluid dynamics which 
respects an SO{d— 1) x 5*0(1, 1) x Z2 symmetry (with d > A) which generalizes the 5*0(3)- 
invariant flow of pLj. The 5*0(1, 1) x Z2 symmetry is the standard boost invariant symmetry 
often considered in the heavy ion literature, and the SO{d — 1) symmetry is a subgroup 
of the 50 (d, 2) conformal symmetry. Our method of constructing the flow involves a Weyl 
rescaling and coordinate transformation of the metric which promoted conformal symmetries 
on R'^"^'^ to isometries on dSd-i x R. In the inviscid case we also found an 5*0(2)-invariant 
flow in three dimensions and a general solution with a broken Z2 symmetry. It should 
be straightforward to use the same techniques to construct flows with different 5*0(1, 1) 
subgroups of S0{d,2) i.e., flows which are not boost invariant. 

When working in the dSd-i x R conformal frame the flow is static in the sense that the 
velocity field has a component only in the de-Sitter time direction. The simplicity of this 
solution allowed us to study linear perturbations around it by decomposing the perturbations 
into scalars and vectors on the 5^ in dS^. In the inviscid case we found analytic expressions 
for the perturbation. The viscous case was studied numerically and in the short wavelength 
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or long time limit. An analysis of the linear perturbations revealed that the 5'0(3)-invariant 
flow suffers from several instabilities associated with a breaking of the Z2 symmetry, i.e., 
instabilities which generate a non-trivial velocity in the rapidity direction. Most of these 
instabilities are cured when the viscosity is small but non-vanishing, and the one that does 
remain occurs at early enough times that it is irrelevant for a study of heavy ion collisions. 
In short, the upshot of our analysis is the the S'O (3)-invariant flow is stable to perturbations 
for parameters similar to the ones reached experimentally at RHIC. 

Finally, following recent interest in higher order moments of the flow of real quark-gluon 
plasmas [31 IH El El El IE] we gave a graphical representation of our analytic expressions for 
the evolution of the second and third moments of the linearized scalar perturbations of the 
flow. In addition, we discussed how our decomposition of the modes in terms of spherical 
harmonics compares with the cumulant expansion of |8j. 
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A de-Sitter space coordinates 



In section we discussed various coordinate systems for de Sitter space. In what follows 
we will provide a slightly more detailed exposition of such coordinate systems emphasizing 
their geometric aspects. Most of this material exists in the literature and can be found in, 
for example, [33]. We have reproduced it for the sake of completeness. In what follows we 
discuss the dS^ geometry most relevant to this work. 

Three dimensional de Sitter space can be visualized as the hyperboloid 

- {Xy + {Xy + {Xy + (X3)2 = (182) 

in R^'^ with line element 

ds^ = -{dXy + {dXy + {dXy + {dXy . (183) 
In what follows we will set L = 1 for convenience. The coordinate system (38), 



X^ = sinh p X^ = cosh p sin 9 cos = cosh p sin 6 sin = cosh p cos 9 

(184) 

covers the hyperboloid. In these coordinates the line element takes the form 

ds"^ = -dp^ + cosh^ p {de^ + sin^ Od^^) . (185) 

The constant p spatial hypersurfaces are two-spheres S"^, parameterized by 6 and (p, whose 
radius shrinks as one goes from p = — ootop = and then expands as p grows from to 



oo. At p = the radius of the S"^ reaches its minimal value, 1. See figure 10 
The coordinate transformation 

2t t 



:i86) 



is identical to (37) except that we have set q = 1 for convenience. Since 

- = X^ - X° , (187) 
r 

spacelike surfaces of constant r are surfaces of constant X^ — X'' on the hyperboloid. Since 
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Figure 10: The dS^ Hyperboloid (with the direction suppressed). Hypersurfaces of 
constant p are two spheres, which are depicted as black circles in the figure. The radius 
of the S"^ is minimized at p = 0. Lines of constant and 6 form the boundaries of the 
hyperbola. The dashed red lines are surfaces of constant r 



r > 0, the (r, x , x ) coordinate system covers only half of the hyperboloid. This is depicted 



in figure 10 



Another way to understand the relation between the (r, a;^,x^) coordinates and the 
{p,9,(j)) coordinates is to note that the spacelike surface of constant p = po corresponds 
to the surface 



= (x'^Y + {x'^Y = (r — (cosh po — sinh Po))(t — (cosh po + sinh po)) 



implying that coshpo + sinhpo < r < oo. Implementing (188) in the relations 



X" , , 1 — r + r 



tan (j) = — cot I 
x^ 



2r 



;i89) 



one finds that 



< < 27r sinh po < cot 9 < oo . 



(190) 



Put differently, at time p = po, the (r, a;^,x^) coordinate system covers the part of the unit 
two- sphere given by 

(x^ + {x'^Y + {x^Y = 1 x^>tanhpo (191) 
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B Conventions 



Our conventions for spherical harmonics are the standard ones implemented in Mathematica 
7.0. The first few are: 



^0,0 

^2,-2 

^2,0 

^2,2 
^3,-3 
^3,-2 

^3,0 

^3,2 
^3,3 



2^ 
2 V 27r 

2 V 27r 

1 fl5 
iV 2^ 
1 fib 

1 /5 
4 

1 /15 



e-"^sin^ 
cos^ 

e^-^sin^ 
e-^^-^sin^^ 
e"'"^ sine cos e 
y?(3cos2e- 1) 



2 V 27r 

1 flh 
iV 2^ 
1 /35 

TT 




e^'^ sine COS e 
e^^-^sin^e 
e-^^-^sin^e 



1 /105 



4 V 27r 
1 /21 



sin 9 cos e 



y^e-^'^sine(5cos^e- 1) 
y^(5 cos^e-3cos e) 
^y^e^^sine(5cos2e-l) 



1 /105 
iV 2^ 
1 /35 



e^'-^sin^ecose 
e^^^ sin=^ e . 



(192) 
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